ABSTRACT 

By  two-phase  flow  theory,  a  new  mathematical  model  of  the  ice  dynamics  in  Arctic  Ocean 
is  formulated.  The  theory  is  developed  from  both  microscopic  and  macroscopic  approaches. 
The  effective  two-dimensional  fundamental  equations  of  ioe  circulation  in  Arctic  Ocean  are  de¬ 
rived  and  compared  with  other  models  used  In  literature.  Our  model  is  (he  most  general  one 
in  literature  and  it  gives  many  new  insights  of  various  important  terms  in  the  compactness 

equation  and  in  the  equations  of  motion  of  ice  dynamics.  Further  development  and  application 
S 

of  thl$  model  to  specific  cases  in  Polar  Ocean  are  recommended.  * 
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I.  Introduction 


For  the  research  of  the  dynamics  of  the  peek  toe  In  the  Arctic  Ocean  and  its  surrounding 
seas,  one  of  the  main  purposes  Is  to  coordinate  theoretical  and  field  investigations  such  that 
die  forecasting  of  ice  motion  (W1  ttmann-MacDonell,  1964)  or  knowledge  of  the  general  circul¬ 
ation  in  the  Arctic  Ocean  and  its  surrounding  seas  (Campbell,  1965),  may  be  improved.  Con¬ 
siderable  theoretical  research  works  have  been  carried  out  (AIDJEX  Bulletins  No.  2  and  No.  3). 
But  our  knowledge  on  the  dynamics  of  pack  ise  is  still  meager  and  more  theoretical  analyses 
are  needed. 

In  the  Arctic  Ocean  and  its  surrounding  seas,  an  enormous  number  of  ice  floes  of  various 
sizes  and  shapes  are  floating  on  the  surface.  Most  of  these  ice  floes,  under  the  influence  of 
winds,  ocean  currents  and  tidal  waves  are  in  constant  but  chaotic  motion  and  are  continuously 
being  deformed.  However,  for  many  practical  reasons,  we  are  Interested  in  the  average  motion 
of  sea  ice  over  a  large  area  such  as  several  tens,  or  hundreds  kilometers  radius  or  larger, 
and  for  an  extended  period  of  time,  say  several  days,  months  or  years.  For  such  large  spatial 
areas  and  long  period,  the  Individual  properties  of  ice  Hoes  become  less  important  and  the 
pack  ice  may  be  considered  as  a  continuum.  In  this  report,  a  new  mathematical  model  for  the 
dynamics  of  pack  ice  in  Arctic  Ocean  and  its  surrounding  seas  is  proposed  so  that  the  pack  ice 
may  be  treated  as  a  continuum. 

The  mathematical  model  of  dynamics  of  pack  ice  may  be  divided  into  two  classes: 

In  the  first  class  which  is  used  by  the  majority  of  the  research  workers  in  this  field,  one 
considers  the  pack  ice  as  a  pseudo-fluid  (Campbell,  1965;  Fel'zenbaum,  1958;  Doronin,  1970), 
while  in  the  second  class,  one  considers  the  pack  ice  as  an  elastic-plastic  material  (Coon  et 
al,  1974).  Moe  t  of  the  computation  of  ice  drift  in  Arctic  Ocean  have  been  carried  out  by  the 
first  class  of  mathematical  model.  Our  model  is  a  most  advanced  mathematical  model  of  the 
first  class. 

Our  main  concept  of  the  mathematical  model  is  the  two-phase  flow  approach.  We  consider 

the  Arctic  pack  ice  as  a  mixture  of  solid  (loe)  and  a  fluid  (sea  water).  Let  the  depth  of  the  ice 

be  H(x,y,  t)  which  is  in  general  a  function  of  x  and  y,  the  two  spatial  coordinates  on  the  surfaoe 

of  the  ocean,  and  the  time  t.  The  maximum  thickness  of  the  loe  H  is  of  the  order  of  ten 

o 

meters.  We  choose  a  typical  length  of  the  ocean  L  to  represent  the  dimension  of  the  surfaoe  of 
the  ocean  which  we  are  Interested  in  suoh  as  x  ~L  and  y~L.  The  representative  length  L  may 
be  of  the  order  of  1,000  kilometers.  Hence  L  is  much  larger  than  Hq.  We  may  divide  the  ice 
packed  ocean  into  two  layers.  In  the  upper  layer,  l.e. ,  z  £  Hq,  where  z  is  the  coordinate  per¬ 
pendicular  to  the  surface  of  the  ocean  measured  from  the  surfaoe  of  the  loe  and  positive 
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downward,  we  have  a  mixture  of  ioe  floes  and  sea  water.  Since  L»H  ,  the  Ice  floes  may  be 

o 

considered  as  small  solid  particles  in  a  fluid  (sea  water).  In  our  mathematical  model,  we  use 
the  concept  of  fluidization  of  solid  particles  in  a  fluid  (Pai,  1971).  Hence  these  ioe  floes  may 
be  considered  as  a  pseudo-fluid.  For  z  -  Ho>  we  have  a  mixture  of  two  fluids;  one  is  the  sea 
water  and  the  other  is  the  pseudo-fluid  of  ioe  floes.  We  may  extend  the  theory  of  two-phase 
flows  of  a  mixture  of  small  solid  particles  and  a  fluid  to  study  the  fundamental  equations  of  ioe 
dynamics  for  the  pack  ice.  We  derive  rigorously  the  fundamental  equations  of  ice  dynamics 
from  the  two-phase  flow  theory  (Pai,  1971),  particularly  the  equation  of  continuity  of  pack  ice 
in  which  the  source  term  depends  on  the  freezing  and  melting  of  ice  floes  and  which  is  coupled 
with  the  energy  equation. 

It  should  be  noted  that  the  three-dimensional  fundamental  equations  just  mentioned  are  not 
the  fundamental  equations  of  ice  circulation  used  by  many  authors  in  the  literature  (Campbell, 
1965;  Fel'  zenbaum,  1958;  Doronin,  1970).  The  fundamental  equations  used  in  literature  by 
these  authors  are  the  integrated  forms  of  the  two-phase  flow  equations  of  three  dimensions. 
Because  of  the  fact  that  L  »  H^,  we  may  Integrate  the  three-dimensional  equations  of  ice  dyna¬ 
mics  with  respect  to  z,  the  vertical  coordinate,  with  the  limits  from  z  =  0  to  z  =  H(x,y,  t).  The 
resultant  equations  depend  only  on  two  spatial  coordinates  x  and  y  and  time,  t,  with  the  thick¬ 
ness  of  ice  H(x,y,  t)  as  an  unknown  parameter  to  be  determined.  This  procedure  is  similar  to 
the  integral  method  which  has  been  extensively  used  in  boundary  layer  flow  (Pai,  1956).  Even 
though  these  integrated  equations  have  been  generally  used,  it  seems  to  us  that  no  one  has  pre¬ 
viously  presented  the  fundamental  equations  of  the  effective  two-dimensional  ioe  flow  in  this 
manner. 

The  integrated  equations  of  motion  of  ice  floes  are  of  a  form  similar  to  the  equations  of 
motion  for  a  single  floe  used  by  Campbell  (Campbell,  1965)  and  others  (Fel'zenbaum,  1958; 
Doronin,  1970;  Rothrock,  1970).  But  the  integrated  equations  give  us  much  more  information. 
In  the  integrated  equations  of  motion,  the  shearing  stresses  on  the  ice  surfaces  z  =  O  and  z  = 

H  occur  such  that  at  z  =  O,  the  shearing  stresses  are  equal  to  the  wind  stresses  Ta  while  at 
z  =  H,  the  shearing  stresses  are  equal  to  foe  water  current  stresses  t  .  Conceptually,  the 
Integrated  forms  of  the  fundamental  equations  are  consistent  with  foe  mathematical  model  of 
pack  ice.  We  cannot  justify  the  derivation  of  foe  equations  of  motion  of  a  continuum  model  of 
pack  ice  by  using  the  equations  of  motion  of  a  single  ice  floe.  From  the  eqxatlons  of  motion 
for  a  single  ice  floe,  we  cannot  obtain  foe  effect  of  compactness.  Furthermore,  we  have  a 
better  understanding  of  various  terms  in  the  Integrated  forms  of  the  fundamental  equations 
such  as  foe  source  terms  and  the  interaction  terms.  We  also  expect  some  new  information 


about  the  effective  two-dimensional  equations  such  as  the  interaction  forces  between  ioe  floes 
and  the  effects  due  to  the  beat  transfer  and  mass  transfer. 

For  z^H,  our  equations  of  two-phase  flow  will  be  reduced  to  those  of  ocean  water  where 
the  compactness  is  zero.  Hence  we  need  only  to  derive  the  fundamental  equations  for  the  upper 
layer,  z  =  H.  Those  in  the  lower  layer  where  z  >H  are  a  special  case  of  those  of  the  upper 
layer  with  zero  compactness. 

We  shall  apply  the  modern  technique  of  fluid  mechanics  from  both  the  macroscopic  (the 
continuum  theory)  and  the  microscopic  (the  kinetic  theory)  point  of  view  to  derive  the  funda¬ 
mental  equations  of  two-phase  flow  of  a  mixture  of  ice  and  water.  Special  attention  will  be 
given  to  the  effect  of  heat  transfer  and  mass  transfer  between  ice,  water  and  air. 

From  the  kinetic  theory  of  gases,  the  fundamental  equations  of  fluid  dynamics  may  be  de¬ 
rived  as  the  transfer  equations  of  the  Boltzmann  equation  of  a  gas  and  some  insights  of  various 
terms  by  these  equations  may  be  obtained.  The  continuum  theory  of  fluid  dynamics  cannot 
give  us  any  insights  of  these  terms.  On  the  other  hand,  because  of  the  mathematical  and  phy¬ 
sical  complexity,  we  can  only  deal  successfully  with  Boltzmann  equation  for  simple  kinetic 
picture  and  for  very  simple  configurations.  For  practical  problems,  we  have  to  use  the  mac¬ 
roscopic  variables  and  the  continuum  theory.  We  postulate  the  fundamental  equations  of  fluid- 
particle  system  based  on  the  conservation  laws  of  mass,  momentum  and  energy.  These  fun¬ 
damental  equations  should  be  consistent  with  those  derived  from  the  kinetic  theory.  We  are 
especially  interested  in  those  new  terms  of  interaction  between  the  fluid  and  solid  particles, 
which  may  be  obtained  by  the  kinetic  theory  but  may  be  ignored  in  the  ordinary  derivation  of 
the  fundamental  equations  by  the  continuum  theory.  Since  the  complete  kinetic  theory  of  the 
fluid-particle  system  has  not  yet  been  developed,  it  is  appropriate  to  use  both  the  macroscopic 
and  microscopic  approaches  to  derive  a  general  set  of  the  three-dimensional  fundamental 
equations  of  ice  dynamics.  From  these  three  dimensional  equations,  the  effective  two-dimen¬ 
sional  equations  of  pack  ice  may  be  obtained. 

With  these  effective  two-dimensional  equations  of  pack  ice,  we  calculate  the  ioe  circula¬ 
tion  in  the  polar  Ocean  and  compare  these  results  with  those  in  the  literature  and  with  the 
field  observations.  We  also  compare  our  model  with  AEDJEX's  model  for  the  similarities  and 
the  differences,  and  briefly  discuss  how  to  use  AIDJEX's  experimental  data  to  verify  the 
essential  points  of  our  model. 

n.  Microsoopic  Approach:  Klnetto  Theory  of  Ioe  Dynamios. 

Since  the  kinetic  theory  of  gases  gives  us  die  Navier-Sbokes  equations  of  a  viscous  fluid 
as  a  first  approximation,  we  would  expect  that  the  klnetto  theory  of  a  gas-particle  system 
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(Pal,  1971;  Kuentzmann,  1973)  would  give  the  fundamental  equations  of  loe  dynamics  (a  mix¬ 
ture  of  ice  floes  and  sea  water)  as  a  first  approximation.  We  consider  an  ideal  case  of  a  gas 
(which  may  represent  the  water  in  our  final  result)  and  many  small  solid  particles  (which  may 
represent  the  elementary  ice  particles  as  we  will  define  later).  The  kinetic  theory  of  this  gas- 
particle  system  may  be  treated  by  using  a  single  particle  distribution  function  for  the  gas, 
fg,  and  that  for  the  solid  particles,  fp.  Each  of  these  distribution  functions  is  governed  by  a 
Boltzmann  equation. 

The  distribution  function  for  the  solid  particles,  fp,  is  different  from  that  of  gas  mole¬ 
cules,  fg,  because  we  have  to  consider  the  different  sizes,  shapes  and  physical  properties  of 
the  solid  particles.  In  general,  the  shapes  of  the  solid  particles  are  arbitrary.  However,  if 
we  were  to  consider  the  arbitrary  shapes  of  the  solid  particles,  the  results  would  be  too  com¬ 
plicated  and  too  detailed.  Hence  we  will  assume  that  the  shapes  of  these  solid  particles  are 
similar.  Thus,  we  need  to  use  one  of  its  dimensions  as  its  size.  With  these  approximations, 
the  distribution  function  of  the  solid  particles  may  be  defined  as  follows: 

The  number  density  of  solid  particles,  dnp,  with  the  size  range  h^  and  h{  +  dh.,  in  the  vol¬ 
ume  x  and  x  +  dx  ,  having  the  particle  instantaneous  velocity,  c  ,  in  the  range  of  c  and  c  + 

J  J  J  J  J  J 

dc.;  and  the  instantaneous  temperature  6  of  the  particle  in  the  range  of  6  and  ©  +  d©  is 

dnp  =  fp(x.,  c^  ©,  h},  t)  d3Cj  d©  dt^  (1) 

where  np  is  the  number  density  of  the  solid  particles,  x^  is  the  jth  spatial  coordinate  of  the 

point  considered,  c  is  the  jth  component  of  the  instantaneous  velocity  of  the  solid  particle  in 
J  3 

the  direction  of  x^,  d  c^  is  the  elementary  volume  of  the  velocity  space,  ©  is  the  instantaneous 
temperature  of  a  solid  particle  which  may  be  varied  from  particle  to  particle,  and  h(  is  the 
length  scale  of  the  particle  which  represents  the  size  of  the  particle  and  which  may  be  different 
for  different  particles.  For  ice  dynamics,  we  should  choose  the  scale  of  such  that  an  ele¬ 
mentary  ice  particle  would  have  uniform  physical  properties  such  as  salinity,  porosity,  etc. 
Hence  the  elementary  ice  particle  in  the  microscopic  analysis  is  much  smaller  than  an  ice 
floe.  In  fact,  an  ice  floe  may  be  considered  as  a  group  of  a  large  number  of  our  elementary 
solid  ice  particles.  The  rate  of  change  of  the  size  of  our  elementary  ice  particle  depends  on 
the  physical  and  mechanical  processes  such  as  evaporation,  condensation,  freezing,  melting, 
disintegration,  hummocktng,  etc.  Hence  is  a  random  variable  in  the  microscopic  analysis. 

The  distribution  Amotion  of  the  gas,  f  ,  may  be  considered  as  a  special  case  of  f  ,  in 

8  P 

which  the  variation  of  die  distribution  function  with  size  parameter  h^,  and  with  temperature  © 


may  be  dropped.  Henoe  from  now  on,  we  shall  consider  only  the  distribution  function,  fp, 
except  where  specified  otherwise. 

The  distribution  function,  f  ,  is  governed  by  the  Boltzmann  equation  which  shows  that  the 

P 

total  rate  of  change  of  f  with  time  must  equal  the  change  of  die  number  of  particles  per  unit 

P 

volume  by  collision  effects  in  the  range  of  variation  considered,  i.e. , 


af  Be  f  ,  90,f  .  /Q.f  \  .  >M  f  v  /if  \  /if  \ 

_IL  +_iiL  J_B.+  _a_/_h£ )+_!-/  PJP.WjL)  +(__£.J  (2) 

at  ax,  m  ac.  aeVmc/  ahVkm  J  \  it  /  \it  /  g 

j  p  j  p  S  i'  P'  p 

where  0^  is  the  jth  component  of  the  body  force  on  a  solid  particle  of  size  ty  and  mass  mp, 
which  is  a  function  of  h{;  is  the  heat  transfer  rate  from  a  particle  of  size  ty,  cg  is  the 
specific  heat  of  the  particle;  is  the  mass  transfer  rate  of  the  particle  of  size  h^  which  may 
be  due  to  the  thermal  and  mechanical  effects,  and  k  is  a  constant  depending  on  the  shape  of  the 
particle.  If  the  particle  is  of  spherical  shape  and  is  its  radius,  k  =  3. 

The  right-hand  side  terms  of  (2)  are  the  collision  terms.  The  first  term  is  the  collision 
between  solid  particles  and  the  second  term  is  the  collision  between  gas  molecules  and  a  solid 
particle  of  the  size  iy  We  shall  study  these  two  terms  in  detail  because  they  have  not  been 
previously  examined  thoroughly.  We  are  especially  Interested  in  the  study  of  collisions  be¬ 
tween  ice  floes  such  as  hummocklng  and  rafting.  The  first  collision  term  between  solid  par¬ 
ticles  has  been  briefly  studied  by  Kuentzmann  (1973).  But  for  ice  dynamics,  many  additional 
effects  such  as  hummocking,  etc. ,  should  be  Included. 

The  collision  terms  for  gas  molecules  and  solid  particle  may  be  split  into  two  parts:  one 
is  due  to  the  mean  flow  of  the  gas  molecule,  l.e. ,  the  fluid  is  considered  as  a  continuum,  and 
the  other  is  due  to  the  random  motion  of  the  gas  molecules.  Thus 


(-B-)  =(^~)  +(-&) 

\  it  /  g  \{t/gm  \5t/j 


the  first  term  (if  /St)  ,  may  be  expressed  as  a  body  force  on  die  solid  particle,  0..  Hence 
p  gm  J 

we  should  use  the  continuum  concept  to  find  out  the  proper  body  force  which  contain  a  drag 
force  similar  to  the  Stokes  formula,  but,  also  including  other  effects  such  as  volume  fraction, 
Reynolds  number,  pressure  gradient,  and  other  effects.  (Pal-Hsleh,  1973V  This  is  the  area 
where  simultaneous  studies  of  the  microscopic  and  die  macroscopic  point  of  view  would  yield 
useful  information. 

The  other  use  of  (2)  is  to  derive  die  fundamental  equations  for  the  macroscopic  approach 

from  the  transfer  equations.  Before  we  derive  the  macroscopic  fundamental  equations,  we 

must  define  the  macroscopic  variables  in  terms  of  the  distribution  function,  f  ,  as  follows: 

5  P 


The  macroscopic  variables  are  the  moments  of  the  distribution  function.  We  have 
(i)  For  zeroth  moment,  we  obtain  the  number  density  and  the  mass  density. 


The  number  density,  n^,  of  the  solid  particles  is 

n  -((ft  (x.,  c  ,  0,  h  ,  t)  d3c  d6  dh 
P  JJJ  P  J  J  1  J  i 

and  the  mass  density,  Pp .  of  the  pseudo-fluid  of  solid  particles  is 


(4) 


p  = 

P 


=ff(n  (h  )  f  (x.,  c  ,  ©,  h,,  t)  d3c.  d0  dh 

JJJ  P  i  p  j  j  *  J  * 


=  m  n  =  p  Z 
P  P  sp 


(5) 


where  P^  =  PgpZ  is  the  partial  density  of  the  pseudo-fluid  of  the  solid  particles  in  the  mix¬ 
ture,  m^h^  Is  the  mass  of  a  particle  of  size  h^,  is  the  mean  mass  of  solid  particles  in 

the  pseudo-fluid,  p  is  the  species  density  of  the  particle,  i.e. ,  the  density  of  the  particle  it- 
sp 

self  without  considering  it  as  a  part  of  the  mixture,  and  Z  is  the  volume  fraction  of  the  solid 
particles  in  the  mixture.  In  two-phase  flow,  the  partial  density  should  be  used  as  we  shall 
show  later. 

(il)  For  the  first  moment  of  the  distribution  function,  we  have  the  following  macroscopic 
variables; 


The  temperature  of  the  solid  particle,  T  ,  is 


T  = 

P 


fm  (h  )  0  f  d  c  d©  dh 


(6) 


P  P 


The  jth  component  of  the  flow  velocity  of  the  solid  particles  is 


V 


"V  °i 

P  P 


f  d  c,  d©  dh, 
P  J  * 


(7) 


(iii)  For  the  second  moment,  we  have  the  jkth  component  of  the  stress  tensor  of  the  pseudo¬ 
fluid  of  the  solid  particles: 


sPjk  -fff%  <hi)  <oj  -  v  <ck  ■  vk»  fP  “S ae  dhi 


(8) 


The  stress  tensor  Sp^  consists  of  the  partial  pressure  p^  of  the  pseudo-fluid  of  solid  particles 
and  the  viscous  stress  tensor  of  the  pseudo-fluid  of  solid  particles  due  to  the  random  montlon 
such  that 


S  =  p  5  -  r 

pjk  p  jk  pjk 


(9) 


where 


=  1  if  j  =  k  and  6 
jk  jk 


=  O  if  j  A  k. 


The  jkth  component  of  the  viscous  stress  tensor  is 


Tpjk* 

Similarly,  we  may  define  other  macroscopio  variables  such  as  the  hest  flux  due  to  con¬ 
duction,  the  interaction  forces  exerted  upon  solid  particles  by  the  pseudo-fluid,  work  done  on 
the  fluid  by  the  entire  particle  cloud,  etc. ,  as  shown  in  reference  by  Pai  (1973). 

The  equation  of  continuity  of  a  species  in  the  mixture  may  be  obtained  by  taking  the  zeroth 

moment  of  the  Boltzmann  equation.  We  multiply  (2)  by  m  (h  )  and  integrate  the  resultant 

3  p  1 

equation  with  respect  to  d  c.  over  the  whole  velocity  space,  with  respect  to  ©  over  the  whole 
temperature  space,  and  with  respect  to  h^  over  the  whole  size  space.  We  may  write  the  re¬ 
sultant  equation  of  continuity  of  the  pseudo-fluid  of  solid  particles  (pack  ice)  as  follows: 


dP 


_i  + 

dt 


*?iui  + 


»vi  + 


a?iwi  = 


O  i 


dx 


dy 


dz 


(10) 


where  p^  =  p^  =  m^  =  P  Z{,  is  the  species  density  of  the  pack  ice,  and  p.  is  the  partial 

density  of  the  pack  ice  in  the  mixture  of  ice  and  sea  water.  We  have  two  different  definitions 

of  density  for  each  species  in  the  mixture.  Let  V  =  V,  +  V  be  the  elementary  volume  of  the 

i  w  ' 

mixture  where  is  the  volume  occupied  by  the  ice,  and  is  the  volume  occupied  by  the  sea 

water.  Let  M  =  M  +  M  be  the  total  mass  of  the  mixture  of  ice  and  water  in  V  where  M  is  the 
i  w  i 

mass  of  ice  in  V  and  is  the  mass  of  water  in  V.  The  species  density  of  ice  is 


M. 


=  °i  <  V  v  V  c > 


(11) 


where  p ^  is  in  general  a  function  of  ioe  temperature  T^,  ice  pressure  p^,  ice  salinity  St,  and 
ice  porosity  £.  The  partial  density  of  the  ice,  p^,  in  the  mixture,  which  is  the  value  of  den¬ 


sity  used  in  the  two-phase  flow  theory,  e.  g. ,  in  (10),  is 


p,  - 


M 


i  = 


M 


l 


V 

_ i 

V 


=  Plzl 


(12) 


where 


V 


l  =  volume  fraction  of  ice  in  the  mixture 


(13) 


The  source  term,  a  ,  depends  on  the  mass  transfer  M  and  the  integral  of  the  collision 
*  P 

terms.  We  shall  study  this  berm  in  detail  in  the  near  future. 

The  equation  of  motion  is  associated  with  the  first  moment  of  the  Boltzmann  equation.  We 

Q 

multiply  (2)  by  m^  and  integrate  the  resultant  equation  with  respect  to  d  c^  d©  dh(  and  obtain 
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the  equation  of  motion  of  the  pseudo-fluid  of  the  pack  ice  as  follows: 


P|ui  +*Wi  +!ViVl  *Wi  -  zip(tvi  + 


IT  a*r 

+  2Z.  p  w.  ft  cos©  =  -  z.  +  ixy  +  i 
1  1  1  1  dx  dx  a: 


+  5Tlxz  + 


+  V^K^)^-^)  (14) 

where  f  =  2  ft  sin©  =  coriolis  force  factor,  ©  is  the  latitude,  and  ft  is  the  angular  velocity  of 
the  Earth. 

The  factor  K(ZJ  is  the  friction  coefficient  of  ice  floes  in  water  which  is  a  function  of  the 
compactness  or  volume  fraction  Z  .  We  should  study  it  in  detail.  For  a  first  approximation, 
we  may  write 


K(Zi)  =  CZ.  (15) 

where  C  is  a  constant. 

The  factor  F  .  is  the  momentum  transfer  per  unit  mass  due  to  the  source  term,  which  may 
ui 

be  obtained  by  studying  the  mass  transfer  term  and  collision  terms  in  the  Boltzmann  equation. 
For  a  first  approximation,  we  may  assume  that  Fu{  =  u^,  i.  e. ,  new  ice  has  the  same  velocity 
as  the  local  flow  velocity  of  die  pseudo-fluid  of  pack  ice. 

In  a  similar  form,  we  have  the  equations  of  motion  in  the  y-  and  z-directions. 

The  stress  tensor  of  the  pack  ice  with  component  Tlxx*  etc.,  [as  defined  in  (9)  and  (10)] 
is  due  to  the  random  motion  of  the  ice  particles.  Since  the  determination  of  these  stress  tensors 
is  very  difficult  or  nearly  impossible  from  Boltzmann  equation,  particularly  for  the  turbulent 
flow,  we  have  to  postulate  the  expression  of  such  stress  in  terms  of  macroscopic  variables 
such  as  velocity  gradient,  etc.  For  Instance,  Doronln  (1970)  used  the  form  similar  to  the  con¬ 
ventional  expression  of  turbulent  stress.  We  may  also  assume  that  the  pseudo-fluid  of  pack 
ice  is  a  non-Newtonian  fluid  or  even  as  an  elastic-plastic  material  as  AIDJEX  modeling  group 
did  (Coon  et  al,  1974).  The  final  decision  should  be  made  by  checking  with  experimental  re¬ 
sults  as  those  will  be  obtained  by  AIDJEX  experiment.  We  may  use  various  expressions  for 
the  stress  tensor  to  calculate  the  ice  circulation  in  the  Arctic  Ocean  to  find  out  the  difference 


due  to  these  expressions  and  check  them  with  field  observations. 

3 

If  we  multiply  (2)  by  m  c  ©  and  Integrate  over  the  spaces  d  c  d©  dh  ,  we  obtain  the  tem- 

p  s  J  t 


perature  or  energy  equation  of  the  pseudo-fluid  of  solid  particles.  This  temperature  equation 


with  the  kinetic  temperature  equation  of  the  particles  gives  us  the  energy  equation  (Pal,  1973) 
of  the  pack  ice.  We  will  study  this  energy  equation  In  detail  in  the  near  future. 


m.  Macroscopic  Approach:  Effective  Two-Dimensional  Equation  of  Ice  Dynamics. 

In  the  macroscopic  approach,  we  consider  the  mixture  of  solid  particles  and  a  fluid  as  a 

mixture  of  two  fluids:  one  is  the  pseudo-fluid  of  solid  particles  and  the  other  is  the  real 

fluid,  gas  or  liquid.  For  each  species  r  In  the  mixture,  we  would  like  to  know  its  velocity 

vector ~q  with  components,  u  ,  v  and  w  ,  its  temperature  T  ,  its  partial  pressure  p  and  its 
r  r  r  r  r  r 

partial  density  B  •  Hence  for  each  species,  we  have  six  fundamental  equations,  i.  e. ,  one 
r 

equation  of  state,  one  equation  of  continuity,  three  equations  of  motion  and  one  equation  of 
energy.  These  fundamental  equations  may  be  derived  from  thermodynamics  (for  equation  of 
state),  and  conservation  laws  of  mass,  momentum  and  energy  (for  the  other  five  equations). 
These  fundamental  equations  should  be  consistent  with  those  derived  from  the  microscopic 
approach,  i.e. ,  the  various  transfer  equations  of  the  Boltzmann  equation  (2).  This  is  the  point 
where  we  may  compare  the  results  of  the  microscopic  and  the  macroscopic  approaches. 

We  have  shown  that  for  ice  movement  in  the  polar  region,  because  two  dimensions  along 
the  surface  of  the  ocean  are  much  larger  than  the  thickness  of  the  ice,  it  is  convenient  to 
treat  the  ice  movement  as  an  effective  two-dimensional  problem.  As  shown  in  reference  by 
Pai-Li  (1971),  the  correct  way  to  derive  the  fundamental  equations  of  these  effective  two- 
dimensional  equations  is  to  Integrate  the  three-dimensional  equations  such  as  (11)  and  (14) 
etc. ,  with  respect  to  the  vertical  coordinate  z  from  z  =  O  to  z  =  H(x,y,  t),  where  H  is  the 
average  thickness  of  the  pack  ice  over  an  elementary  area  at  the  point  (x,y,t). 

We  are  going  to  show  some  details  for  such  derivations  which  will  give  us  new  information 
for  ice  dynamics. 

(a)  Equation  of  continuity  of  pack  ice  -  compactness  of  pack  ice,  N. 

Equation  (10)  may  be  written  as  follows: 


*Vi  +  aViui  +  aziVi 


+  *Viwi  = 


l 


(10a) 


at  ax  ay  az 

Now  we  integrate  (10a)  with  respect  to  z  from  z  =  O  to  z  =  H(x,y,t)  and  have  the  following 
equation: 


H  H  H  H  H 

f8V i  dz+(8Wl  dz  +  f'Vl'i  ,  +  dz  =  fff,dz  (16) 

j—  y —  j— ■ -  J—  J 1 

Since  the  ice  remains  on  the  top  layer  of  the  ocean,  we  may  assume  that  at  z  =  O  and  z  -  H, 

Wj  =  O  or,  if  we  consider  freezing  and  melting  ice,  the  rate  of  Increase  or  decrease  at  the 
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surface  z  =  O  and  z  =  H  are  approximately  equal,  then  the  fourth  term  on  the  left-hand  side 


of  (16)  vanishes.  We  have  also  the  following  relation 

W  F  «Zl/Dl 


— ~  fz  p.  dz  =  Z  P  + 

0 1 J  i^i  i^i  0t 


dz 


(17) 


at 


and  similar  relations  for  the  second  and  the  third  terms  on  the  left-hand  side  of  (16).  Further¬ 
more,  since  H(x,y,  t)  is  small  with  respect  to  the  horizontal  dimension  L,  we  may  assume 


that  the  variations  of  p^,  Z  ,  u{  and  v^  with  respect  to  z  are  negligibly  small  and  we  may  re¬ 


place  their  values  by  the  average  values  of  these  variables  between  the  limits  z  =  O  to  z  =  H. 


Hence,  we  have  the  following  relations: 
H 


Jv. 


dz  =  H  N  D 


i 


(18a) 


H 


l 


Zi/0iUi  ^  =  H  N/0i  Ui 


1 


ZlPiVi  dz  =H  N/oi  vt 


(18b) 

(18c) 


H 


/- 


dz=  Hff[  =  Ei 


(18d) 


where  N  is  the  compactness  of  the  pack  ice  which  is  the  average  volume  fraction  of  the  pack 
ice  between  z  =  O  to  z  =  H(x,y,  t)  at  any  point  on  the  surface  of  the  ocean  (x,y)  and  at  a  given 
time  t. 

Substituting  (17)  and  (18)  into  (16),  we  have  the  equation  of  continuity  of  pack  ice  in  the 
effective  two-dimensional  flow  of  ice  dynamics  in  polar  ocean  as  follows: 


8NHP1  .  aSHP|U|  +  aNHPlVl  -NP  JVL+e. 
at  ax  ay  1  Dt 


(19) 


where 

D 


_ i 

Dt 


+  u 


i 


8 

8x 


+  v, 


i 


(20) 


8t  ax  ay 

Equation  (19)  may  be  written  as  an  equation  for  compactness  N  as  follows: 


PH  D1N 
1  Dt 


=  -  PjNH( 


aUl  +  aVl  -  NH  DIPi  +  E 


l 


(21) 


0x  ay  Dt 

It  is  interesting  to  compare  our  equation  of  compactness  N,  (21)  with  those  equations  of  com¬ 
pactness  usually  found  in  literature. 

Rothrock  (1970)  gave  the  following  equation  for  compactness: 


P  N  =  -  P  NH(  *Ui  +  Hi  )  -  P.NH*  +  0 

Dt  1  8x  ay  2 


(22) 
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The  differences  between  (21)  and  (22)  are  trl-fold:  (i)  on  the  left-hand  side  of  these  equations, 

H  in  (21)  replaces  N  on  Rothrock's  equation  (22).  It  is  evident  that  N  in  (22)  must  be  a  print¬ 
ing  error  because  the  dimensions  of  this  term  are  not  die  same  as  those  of  all  (he  other  terms 
in  (22).  (11)  NH(D  p  /Dt)  of  (21)  replaces  NH  <|r  of  (22).  Rothrock  gave  little  information 
about  his  function  <|i  .  We  have  no  problem  in  explaining  our  term,  NHfD^  p^/Dt).  (ill)  There 
are  some  basic  differences  in  concept  between  (21)  and  (22).  In  general,  the  velocities  of  ice 
(u.  and  v.)  are  different  from  those  of  water  (u  and  v  )  in  the  top  layer,  O^z^H.  The 
distinction  between  these  velocities  should  be  noted,  and  in  (21),  the  use  of  (he  operator 
(Dj/Dt)  should  be  kept  in  mind.  In  Rothrock's  analysis,  such  a  distinction  was  not  mentioned. 
Rothrock  separated  the  source  term  into  two  parts,  and  in  (22),  only  one  part  of  the  source 
term,  i.e. ,  the  rate  of  formation  of  ice  by  freezing  open  water  is  Included.  In  Rothrock's 
analysis,  if  one  makes  the  distinction  between  the  freezing  at  the  surface  z  =  O  and  freezing 
at  O  <  z  ?  H,  one  cannot  use  the  effective  two-dimensional  approach.  For  an  effective  two- 
dimensional  approach,  the  source  term  should  be  the  average  source  due  to  freezing  or  melt¬ 
ing  at  the  point  (x,y).  As  a  result,  our  E  =  Ha.  is  different  from  Rothrock's  <&  . 

11  ^ 

We  see  that  the  thickness  H  occurs  in  every  term  of  (21)  if  we  write  =  H  o  ,  and  the 
compactness  equation  is  Independent  of  the  thickness  of  fee  ioe.  This  result  is  essential  so 
that  the  pack  ice  may  be  considered  as  a  two-dimensional  continuum. 

Various  Russian  authors  have  used  simplified  equations  for  compactness.  We  are  going 
to  compare  their  simplified  equations  wife  our  exact  equation  of  compactness,  (21). 

Drogaitsev  (1956)  considered  fee  conservation  of  fee  mass  of  ice  per  unit  area  and  obtained 
the  following  equation  for  the  two-dimensional  ice  dynamics. 

aNPlH  +  >1WPlUl  .  aNHVl  >  o  P3, 

at  ax  ay 

Comparing  (23)  with  our  exact  (19),  we  see  that  Drogaitsev  neglected  the  source  term  E^. 
Furthermore,  he  missed  the  Important  term  due  to  fee  variation  of  fee  thickness  of  ioe,  H. 

This  is  fee  main  difference  between  the  effective  two-dimensional  treatment  of  ice  dynamics 
and  the  ordinary  two-dimensional  approach  obtained  by  simply  dropping  the  variation  with  z- 
coordinate.  As  we  shall  show  later,  if  we  simply  drop  fee  z-  variation  In  fee  three-dimensional 
equations  of  motion,  we  will  never  get  fee  terms  due  to  fee  air  and  water  stresses  on  fee  sur¬ 
face  of  the  ice.  Thus,  such  a  simple  equation  of  continuity  (23)  is  not  consistent  wife  the 
correct  two-dimensional  equations  of  motion  of  toe  dynamics.  Therefore,  (23)  is  oorreot  only 
if  we  neglect  the  source  term  and  if  we  assume  that  fee  thickness  of  fee  ioe,  H,  Is  a  constant. 
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Nikoforov  (1957)  and  Doronln  (1970)  used  the  following  equation  for  compactness: 

8Nu  0Nv 

0N_  =  -  (  1  +  1  )  (24) 

at  ax  ay 

If  we  compare  (24)  with  our  (21),  we  see  that  both  Nikiforov  and  Doronln  neglected  the  source 
term  o{,  and  the  variation  of  the  species  density  of  ioe,  p  .  Thus  (24)  is  valid  only  if  we 
assume  that  there  Is  no  source  term  and  that  the  variation  of  the  species  density  of  ice,  p^, 
is  negligible. 

Similarly,  the  three-dimensional  equation  of  continuity  of  the  water  in  the  mixture  is 
a(i-z.)p  ,  a(i-z.)  p  u  a(i-z,)p  v  a(i-z,)  p  w 

i  rw  +  i  rw  w  +  '  l  rw  w  +  i'  Kw  w  = 


Where  the  subscript  w  refers  to  the  value  of  water  in  (he  mixture,  and  we  assume  °  = 

w 

o^,  i.e. ,  we  neglect  the  evaporation  or  condensation  of  water  vapor  in  the  mixture. 

If  we  integrate  (25)  with  respect  to  z  from  z  =  O  to  z  =  H,  and  use  the  same  arguments  as 
in  the  case  of  ice  and  similar  notations,  we  have  the  equation  of  continuity  of  water  in  the 
effective  two-dimensional  flow  of  ice  dynamics  as  follows: 


P  H  °wN  =  (1- 
W  Dt” 


a-N)  PwH^W  *  !!»)  ♦  (l-N)  ,1-N,  (PwWw,h  m 


where 


w  =  a  +  u 


and  ( P„rw,JH  is  016  value  of  (p  w  )  at  z  =  H  which  may  or  may  not  be  negligible.  The  value 

W  W  Jl  w  w 

of  ( p  W,J  at  z  =  O  is  negligible  in  our  case.  The  difference  between  CD,/ Dt)  and  (D  /Dt)  in 
w  w  i  w 

(21)  and  (25)  should  be  noted.  In  general,  u.  is  not  equal  to  u  and  v.  is  not  equal  to  v  .  The 

1  w  i  w 

source  term  Ej  =  in  (21)  and  (26)  are  the  same.  Equation  (26)  shows  that  the  vertical  com¬ 
ponent  of  the  water  at  the  bottom  of  the  ice,  z  =  H  may  affect  the  compactness  of  the  ice.  This 
is  a  point  which  has  not  been  extensively  investigated. 

(b)  Equations  of  motion  of  pack  ice :  Interaction  forces  between  ice  floes. 


Similar  to  the  analysis  of  the  equation  of  continuity,  we  may  obtain  the  three-dimensional 
equations  of  motion  of  ice  dynamics  from  the  first  moment  of  the  Boltzmann  equation  (2).  The 
complete  three-dimensional  equations  of  motion  Including  heat  and  mass  transfer,  and  the 
general  collision  terms  have  not  been  studied.  Some  simple  models  of  two-phase  flow,  how¬ 
ever,  have  been  worked  out  (Pal,  1973;  Kuentzmann,  1973).  We  shall  study  the  general  three- 
dimensional  equations  of  motion  of  ice  dynamics  in  the  near  future.  We  will  now  derive  the 
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effective  two-dimensional  equations  from  the  simple  three-dimensional  equations  of  motion 
for  the  two-phase  flow  of  ice  dynamics  by  integrating  the  three-dimensional  equations  of  toe 
dynamics  (Pal,  1973;  Kuentzmann,  1973)  with  respect  to  z  from  z  =  O  to  z  -  H.  We  are  es¬ 
pecially  interested  in  the  new  information  obtained  in  this  manner  as  compared  with  the  deri¬ 
vation  of  the  effective  two-dimensional  equations  of  motion  from  the  foroe  balance  of  a  single 
ice  floe. 

The  three-dimensional  equations  of  motion  of  the  pseudo-fluid  of  pack  ice  in  the  mixture 
of  ice  and  sea  water  are  given  in  (14). 

New  we  will  derive  the  equations  of  motion  of  the  effective  two-dimensional  dynamics  of 
pack  ice  by  integrating  (14)  with  respect  to  z  from  z  =  O  to  z  =  H.  With  the  same  arguments  as 
in  die  case  of  derivation  of  the  continuity  equation  (19),  we  have  the  following  effective  two- 
dimensional  equations  of  motion  for  the  pseudo- fluid  of  pack  ice: 


PiUi  -  f  v  =  -  1  9p  +  _1_  (  8Tlxx  +  *Tixv)  -  aTlxzo 
Dt  Pi  ax  Npj  ax  ay  NHpj 


+  T  IxzH 

+  K(N) 

(u  - 
'  w 

V 

NHpt 

NPi 

Vl  +fu 

Dt  1 

-  -  1 

PT 

9d  + 

1 

»y 

NPi 

tx 

+  TiyzH 

+  m 

(V  - 

w 

V 

NHpt 

ay 


NHo 


(27) 


(28) 


The  source  terms  disappear  because  we  assume  F  ,  =  u.,  etc.  The  terms  t  and  t 

ui  i  ixzo  iyzo 

are  respectively  the  corresponding  wind  stresses  on  the  top  of  the  surface  of  ioe  floes,  z  =  O, 

in  the  direction  of  x  and  y.  The  terms  T.  „  and  t  _  are  respectively  the  water  stresses 

ixzH  iyzH 

acting  on  die  lower  surface  of  the  ice,  z  =  H,  in  the  direction  of  x  and  y. 

There  are  two  types  of  interaction  foroes  in  (27)  and  (28):  (1)  one  is  the  friction  force  due 
to  the  difference  in  velocities  between  ioe  and  water,  i.e. , 


KfliL 

N  Pj 


<  \  -  V 


(28) 


where  q^  and  are  respectively  the  velocity  vector  of  water  and  toe  in  the  mixture.  This  is 

the  drag  foroe  of  die  ice  floes  moving  In  water,  which  is  mainly  due  to  the  collision  between 

the  mean  water  foroe  and  the  ice  floes  as  we  mentioned  in  H.  This  foroe  was  first  proposed 

by  Sverdrup  (1928)  but  neglected  in  muoh  of  the  literature  (Campbell,  1965;  Doronin,  1970; 

13 


F 


Fel'zenbaum,  1961).  In  Sverdrup's  original  proposal,  the  water  velocity  q_  (u  ,  v  )  has 

W  w  w 

been  neglected.  In  our  general  mathematical  model,  the  water  velocity  components  (u^,  v  ) 

are  retained,  which  may  or  may  not  be  negligible  to  ice  velocity  components  (u^,  v^).  If  the 

ice  motion  is  mainly  due  to  the  wind  stress,  the  velocity  of  the  ice  (u  ,  and/or  v{)  may  be  much 

larger  than  the  water  velocity  (u  and/or  v  ).  This  is  the  reason  why  Sverdrup  neglected  the 

w  w 

water  velocity  in  his  formula.  On  the  other  hand,  if  we  study  the  tidal  effect  of  ice,  the  water 
velocity  due  to  the  tidal  waves  may  be  much  larger  than  that  part  of  velocity  of  ice  due  to  the 
tidal  force.  Thus,  we  have  to  include  the  water  velocity  in  the  formula  of  the  friction  forces, 
(ii)  The  other  interaction  force  is  the  internal  stresses  of  ice  floes  due  to  their  random  motion. 
One  of  the  main  objectives  of  AIDJEX  is  to  find  the  proper  expression  of  these  internal  stres¬ 
ses.  According  to  Doronin  (1970),  we  may  write 

8Tixx  +  aTixv  =_J_(€  Npt^i  )  +  _9_(ctNpi  ^i)  (30a) 

ax  ay  ax  t  ax  ay  ay 

8Tixy  +  aTiyy  =  a  UtNp  *VJ  )  +  a  (^Np  *_i  )  (30b) 

ax  ay  ax  ax  ay  1  ay 

where  e  is  the  turbulent  exchange  coefficient  to  be  determined  by  experiment  or  estimated  by 

V 

current  knowledge  on  turbulent  flow.  One  of  the  new  features  of  Doronin's  formula  is  that  he 
includes  the  effect  of  compactness.  His  formula  is  consistent  with  the  two-phase  approach  be¬ 
cause  in  two-phase  flow  analysis,  the  partial  density  N  p{  should  be  used  in  all  equations  in¬ 
stead  of  the  species  density  P^.  We  think  that  (30)  may  be  used  in  (27)  and  (28)  as  a  first 
approximation,  but  there  is  a  difference  between  our  equations  and  those  of  Doronin.  to  our 
equations  (27)  and  (28),  the  stress  terms  (30)  are  divided  by  Np{,  while  in  Doronin's  equations 
of  motion,  they  are  divided  by  Pj  only.  When  we  consider  the  mixture  of  ice  and  water,  we 
should  use  the  partial  density  Np^  in  the  equations  of  motion.  Hence,  Doronin's  equations  of 
motion  were  not  correct  by  using  the  species  density  p{  in  the  Inertial  terms. 

Similarly,  we  may  derive  toe  effective  two-dimensional  equations  of  motion  of  water  from 
two-phase  flow,  three-dimensional  equations. 

We  shall  also  study  the  three-dimensional  equations  of  energy  for  both  water  and  ice  from 
the  Boltzmann  equation  and  their  effective  two-dimensional  form  for  a  complete  study  of  ioe 
dynamics.  We  may  also  derive  equations  for  salinity  for  the  two-phase  analysts.  The  equation 
of  salinity  of  water  is  the  diffusion  equation  of  salt  in  water.  The  salinity  of  ice  depends  on 
the  salinity  of  water  from  which  the  Ice  is  formed  and  other  processes  (Zubov,  1945),  such  as 
the  rate  of  ice  formation,  the  state  of  sea  during  the  ice  formation,  and  the  age  of  the  ioe. 
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IV,  Comparison  of  our  mathematical  model  with  other  models. 

Now  we  compare  our  model  with  (hose  other  models  in  the  literature  (Campbell,  1965; 
Doronin,  1970;  Fel'zenbaum,  1961;  Coon,  et  al,  1974).  In  general,  there  are  two  different 
types  of  mathematical  models  for  pack  ice.  One  is  the  fluidized  model  in  which  (he  pack  loe 
is  considered  as  a  pseudo-fluid.  Our  model  is  essentially  a  very  generalized  version  of  this 
model.  Those  models  in  references  by  Campbell  (1965),  Doronin  (1970),  Fel'zenbaum  (1961) 
and  Sverdrup  (1928)  belong  to  this  type.  We  are  going  to  discuss  the  relationship  between  our 
model  and  those  other  fluidized  models  in  detail  in  the  following  sections.  The  other  type  is 
the  elastic-plastic  model  in  which  the  pack  ice  is  considered  as  an  elastic-plastic  material 
known  as  AIDJGX's  model  (Coon,  et  al,  1974).  We  shall  discuss  the  similarities  and  differ¬ 
ences  of  our  model  and  the  AIDJEX*s  model  later. 

a.  Fluidized  Model. 

In  our  two-phase  flow  model,  we  consider  the  mixture  of  loe  and  water  as  a  mixture  of 
two  fluids  that  have  the  following  variables: 


<?  ,  T  ,  p  ,  P  ,  S 
r  r  r  r  r 


(I) 


The  species  r  may  be  the  pseudo-fluid  of  pack  ice  (r  =  1)  or  the  water  (r  =  w);  q  is  the  mean 
flow  velocity  vector  with  component  u,  v  and  w;  hence  there  are  six  velocity  components 
(u  ,  v  ,  w  ,  u  ,  v  ,  v:  ).  There  are  two  temperature  (T  and  T  ),  and  two  pressures  (p 

111  W  W  W  1  W  1 

and  p  )  which  may  be  expressed  in  terms  of  die  total  pressure  of  the  mixture  p  and  the  vol- 

W 

ume  fraction  Z^,  i.e. ,  pt  =  Z^p  and  p^  =  (1  -  ZJ  p.  There  are  two  partial  density  P{  and  p^, 
which  are  respectively  Pj  =  ZjPj  and  Pw  =  (1  -  Z{)  p^.  The  species  density  (pt  or  p^)  is 
given  by  the  equation  of  state,  i.e. , 


Pr  -  pr  <Tr.  Pr.  Sr>  (31) 

where  Sr  is  the  salinity  of  the  species  r. 

In  the  most  general  case,  we  should  consider  the  14  variables  represented  by  (I).  But 

because  of  the  fact  that  the  dimensions  along  the  ocean  surface,  L,  are  much  larger  than  the 

thickness  of  the  ice,  we  may  consider  the  effective  two-dimensional  flow  of  loe  dynamics  as 

shown  in  Section  m.  In  this  effective  two-dimensional  flow,  we  consider  the  total  pressure 

P  as  a  known  function  given  by  the  atmospheric  pressure  above  die  pack  ice  and  the  inclination 

of  the  ocean  and  the  vertical  velocity  component  w  ,  which  is  much  smaller  than  the  horizontal 

r 

velocity  components  u  and/or  v  ,  may  be  eliminated  by  taking  the  mean  values  across  the 
thickness  of  die  ice  layer.  Thus,  we  have  to  deal  only  with  11  variables,  i.e. , 


"i-  v  v  v  v  v  pt*  pw' N'  sr  sw  m 

Our  model  (II)  is  governed  by  11  fundamental  equations  for  the  11  variables  (II).  These  funda¬ 
mental  equations  are: 

(i)  Equation  of  compactness  (21) 

(ii)  Four  equations  of  motion  (27)  and  (28) 

(iii)  Two  equations  of  state  (31) 

(iv)  Two  equations  of  energy  (T^  and  T  ) 
and  (v)  Two  equations  of  salinity  (S^  and  S^) 

The  authors  believe  that  such  a  complete  analysis  has  not  been  studied  yet.  It  is  our  in¬ 
tention  to  study  such  a  complete  model  eventually. 

We  may  simplify  our  complete  model  (II)  by  various  approximations.  Since  the  salinity  of 
ice  depends  on  the  salinity  of  sea  water  from  which  the  ice  is  formed,  the  rate  of  ice  forma¬ 
tion,  the  state  of  sea  during  the  ice  formation  and  the  age  of  ice,  no  simple  relationship  for 
the  determination  of  ice  salinity  has  been  found  in  connection  with  pack  ice  circulation.  Hence, 
in  most  analysis,  we  may  assume  and  as  constants  or  simple  functions  of  Pf  and  Tf. 

Then  we  have  only  nine  variables: 

u.,  v.,  u  ,  v  ,  T,,  T  »  P.t  P  f  N  (III) 

1  r  w  w  i  w  i  w 

We  feel  that  model  (III)  should  give  good  results  for  ice  circulation  in  the  polar  ocean  by  solv¬ 
ing  the  equation  of  compactness,  four  equations  of  motion,  two  equations  of  energy,  and  two 
equations  of  state  together  with  the  equations  for  air  above  the  ice  and  the  equations  of  water 
below  the  ice. 

Most  of  the  analyses  in  the  literature  are  much  simpler  than  our  model  (III).  Our  model 
may  be  reduced  to  these  simplified  models  as  follows: 

(i)  The  simplest  fluid  model. 

The  simplest  fluid  model  considers  only  the  four  variables: 

u.,  v,,  u  ,  v  (IV) 

iiw  w 

In  references  by  Campbell  (1965),  Fel'zenbaum  (1961)  and  Sverdrup  (1928),  this  model  has 

been  used.  In  order  to  reduce  our  model  (ID)  to  the  simplest  model  (IV),  we  have  to  assume 

that  the  temperatures  T  and  the  denslfy  p  are  constant.  We  also  assume  N  =  constant.  For 

r  r 

Instance,  in  references  by  Campbell  (1965),  Fel’zenbaum  (1961)  and  Sverdmp  (1928),  the 
compactness  is  implicitly  assumed  to  be  unify.  Furthermore  in  these  references,  different 
types  of  equations  of  motion  have  been  used.  In  reference  by  Campbell  (1965),  only  the  viscous 


stresses  (with  N  =  1)  of  (30)  are  used  but  not  the  Interaction  force  (29),  while  in  reference  by 
Sverdrup  (1928),  the  interaction  force  (29)  with  q^  =  O  but  not  the  viscous  stresses  (30)  have 
been  used.  In  reference  by  Fel'zenbaum  (1961),  both  interaction  foroe  (29)  and  the  viscous 
stresses  (30)  are  neglected.  In  our  simplest  model,  we  would  like  to  use  both  interaction 
force  (29)  and  the  internal  stresses  (30)  with  N  =  1. 

(ii)  Doronin  Model. 

In  reference  by  Doronin  (1970),  the  following  variables  are  considered: 


N 


(V) 


U,,  V,,  U  ,  V  , 

i  i’  w’  w 

Doronin  added  an  equation  of  compactness  (24)  to  the  equations  of  motion.  Furthermore, 

Doronin  implicitly  assumed  in  his  model  that  the  velocity  components  of  ice  and  water  in  the 

laver  O  =  z  -  H  are  equal,  i.e. ,  u  =  u  ,  v  =  v  so  that  the  equations  of  motion  of  ice  (27) 

i  w  l  w 

and  (28)  are  solved  with  the  equation  of  compactness  (24).  Because  of  the  assumption  of  u  =  u 

1  w 

and  v  =  v  ,  the  friction  force  (29)  is  identically  equal  to  zero.  Doronin  also  considered  the 
i  w 

thermal  effects  of  compactness  and  the  thickness  of  ice  in  an  empirical  manner  by  assuming 
that  the  temperature  is  known  and  given  by  the  temperature  of  the  atmosphere  at  a  certain 
height.  After  he  obtained  the  dynamic  effects  from  (24),  (27)  and  (28),  he  added  the  thermal 
effects  to  his  results  by  some  empirical  relationships. 

We  may  improve  Doronin' s  model  (V)  by  considering  the  following  variables: 


u.,  v.  u  ,  v  ,  T,  N 
1  l  w  w  ’ 


(VI) 


where  T  =  which  may  be  different  from  the  atmospheric  temperature  and  which  is 

given  by  the  energy  equation  of  the  mixture.  Furthermore,  we  may  calculate  the  case  u  A  u 

1 

and  v.  /  v  in  the  top  layer  O  =  z  =  H. 

I  w 

b.  Comparison  of  Pai-Li  Model  with  AIDJE^s  Model. 

The  AIDJEX*s  dynamic  model  and  its  relationship  to  the  atmospheric  and  ocean  global 
models  are  shown  in  Figure  1,  Figure  6  of  reference  by  Untersteiner  (1974).  Our  mathematical 
model  is  shown  in  Figure  2.  There  are  some  similarities  and  some  differences  between  these 
two  models. 

The  similarities  of  these  two  models  are  as  follows: 

In  both  models,  the  essential  efforts  are  given  to  the  fondamental  equations  which  govern 
the  properties  of  the  pack  ice  as  a  continuum,  such  as  (he  ice  flow  velocity,  ioe  stresses,  ice 
thickness,  etc.  The  fondamental  equations  of  the  ioe  layer  should  be  solved  with  the  flutd 
dynamics  of  the  atmosphere  above  the  ioe  and  that  of  the  ooean  current  below  the  ice.  Thus, 
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Atmosphere/ocean 


Figure  1.  The  AIDJEX  dynamic  ice  model  and  its  relationship  to  the  atmosphere  and  the 
ocean  global  models.  The  nature  of  the  expected  model  output  and  potential 
applications  to  practical  problems  are  also  Indicated. 


APPLICATIONS 
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Figure  2.  Pai-Li  mathematical  model  of  pack  ice  and  its  relationship  to  the 
atmosphere  and  the  ocean  global  models 


In  Figures  1  and  2,  only  the  dynamic  model  part  Is  different. 

The  primary  differences  of  these  two  models  are  as  follows: 

(i)  The  distribution  functions  in  these  two  models  are  of  entirely  different  concept.  In  our 
model,  the  distribution  function  is  connected  with  the  microscopic  analysis  while  in  AIDJEX's 
model,  die  distribution  function  is  connected  with  the  macroscopic  approach.  We  use  the 
Boltzmann  equation  to  determine  the  distribution  function.  In  AIDJEX’s  model,  the  equation  of 
ice  distribution  is  solved  simultaneously  with  the  macroscopic  equation  of  ice  flow  velocity. 

(ii)  In  die  Pai-Li  model,  the  pack  ice  is  considered  as  a  mixture  of  ice  and  water,  and  the 
compactness  of  pack  ice  is  considered  as  a  dependent  variable.  In  the  AIDJEX  model,  the 
pack  ice  is  considered  as  an  elastic-plastic  material.  The  compactness  is  not  considered  as 
one  of  the  dependent  variables. 

(iii)  In  the  AIDJEX  model,  the  constitutive  equation  is  assumed  according  to  the  elastic- 
plastic  material.  In  Pai-Li' s  model,  the  consltitutive  equation  is  that  of  a  viscous  or  turbu¬ 
lent  fluid.  For  instance,  a  constitutive  equation  similar  to  that  of  Doronin  (30)  may  be  used. 
Furthermore,  in  the  Pai-Li  model,  no  limitation  on  the  constitutive  equation  is  imposed.  If 
we  assume  that  the  pack  ice  may  be  considered  as  a  non-Newtonian  fluid,  a  rather  general  con- 
situtive  equation  including  that  of  an  elastic-plastic  material  may  be  used  to  calculate  the  ice 
circulation  in  the  Polar  Ocean. 

(iv)  In  the  AIDJEX  model,  the  thermodynamic  effect  is  Included  by  assuming  that  the 
temperature  of  ice  is  known,  hi  Pai-Li' s  model,  the  temperature  of  the  pack  ice  will  be  cal¬ 
culated  from  the  energy  equation  of  the  pack  toe  using  the  atmospheric  and  oceanic  data  as 
boundary  conditions. 

Since  there  are  several  significant  differences  in  these  two  models,  numerical  solutions 
using  both  models  should  be  carried  out  and  compared  with  each  other  and  with  observational 
data. 

V.  Summary  and  Conclusions. 

The  following  are  the  summary  and  conclusions  of  our  mathematical  model: 

(1)  Based  on  two-phase  flow  theory  (a  mixture  of  a  fluid  and  small  solid  particles),  a  new 
mathematical  model  of  the  ice  dynamics  in  the  Artie  Ocean  is  formulated. 

(2)  This  mathematical  model  is  developed  from  both  the  microscopic  and  the  macroscopic 
approaches. 

(3)  In  the  microscopic  approach,  the  distribution  functions  for  the  pack  toe  and  for  sea 
water  are  studied  and  so  are  the  Boltzmann  equations  which  govern  these  distribution  functions. 
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The  relations  between  the  distribution  Auction  and  the  macroscopic  variables  are  given.  The 
derivations  of  the  macroscopic  equations  from  the  transfer  equations  of  Boltzmann  equation 
are  discussed. 

(4)  In  the  macroscopic  approach,  the  derivations  of  the  effective  two-dimensional 
equations  of  pack  ice  are  discussed  in  detail,  particularly  the  equation  of  compactness,  the 
equations  of  motion  and  the  interaction  forces  between  ice  and  water. 

(5)  We  compare  our  mathematical  model  with  various  models  used  in  literature  and  show 
that  our  model  is  the  most  general  one  for  the  fluidized  model. 

(6)  We  also  show  that  under  various  simplified  approximations,  our  model  may  be  reduced 
to  those  models  in  literature.  We  also  discuss  the  improvement  of  those  simplified  models  in 
literature  based  on  our  mathematical  model. 

VI.  Recommendations . 

The  following  problems  are  suggested  for  further  development  and  application  of  our 
mathematical  model. 

(1)  We  will  apply  some  of  our  simplified  models  such  as  (V)  to  reinvestigate  the  ice  cir¬ 
culation  in  the  central  portion  of  Arctic  Basin.  The  main  point  is  to  show  the  improvements 
of  our  model  (V)  over  those  by  Campbell  (1965)  and  Fel'zenbaum  (1961)  where  the  effect  of 
compactness  is  not  Included.  We  may  also  show  the  effects  of  various  interaction  forces  such 
as  given  by  (29)  and  (30). 

(2)  We  will  apply  our  model  (V)  to  study  the  ice  movement  near  the  coastal  region  or  in 
the  marginal  zone  so  that  the  effect  of  tidal  wave  will  be  included.  This  problem  is  very  im¬ 
portant  from  practical  point  of  view  but  the  numerical  calculation  has  not  been  carried  out. 

(3)  We  shall  study  the  energy  equation  of  the  effective  two-dimensional  ice  dynamics  in 
detail.  After  we  obtain  the  energy  equation,  we  shall  apply  our  model  (VI)  to  study  the  ice 
circulation  in  the  central  portion  of  Arctic  Ocean  and  compare  with  Doronin's  results  (1970) 
with  ours. 

(4)  We  will  apply  our  model  (VI)  to  the  coastal  regions.  We  may  contribute  to  the 
Improvement  of  long  range  ice  prediction  by  considering  the  thermal  effect  and  the  dynamic 
effect  simultaneously. 

(5)  From  the  energy  equation  of  two-phase  flow,  we  will  study  unsteady  toe  dynamic 
problems  such  as  the  growth  and  decay  of  the  young  toe  and  one-year  ice,  particularly 
those  near  the  coastal  region  which  is  a  very  Important  practical  problem. 
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